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ABSTRACT 

Using the 2MASS 2'^'^ Incremental Data Release and the Zodiacal-Subtracted 
Mission Average maps of COBE/DIRBE, we estimate the cosmic background 
in the J (1.25 //m) and K (2.2 fim) bands using selected areas represent- 
ing ~ 550 deg^ of sky. We find a J background of 22.9 ± 7.0 kJy sr^^ 
(54.0 ± 16.8 nW m"^ sr"^) and a K background of 20.4 ± 4.9 kJy sr~^ 
(27.8 ± 6.7 nW m~^ sr~^). This large scale study shows that the main un- 
certainty comes from the residual zodiacal emission. The cosmic background we 
obtain is significantly higher than integrated galaxy counts (3.6 ± 0.8 kJy sr~^ 
and 5.3 ± 1.2 kJy sr~^ for J and K, respectively), suggesting either an increase of 
the galaxy luminosity function for magnitudes fainter than 30 or the existence of 
another contribution to the cosmic background from primeval stars, black holes, 
or relic particle decay. 

Subject headings: cosmology: observations — diffuse radiation — infrared: gen- 
eral 
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1. Introduction 

The cosmic infrared background (CIB) is an important cosmological constraint on the 
star formation history in the Universe (Dwek et al. 1998; Gispert et al. 2000). The first 
detection of the CIB was reported by Puget et al. (1996) for sub-millimetric wavelengths 
using data from the Far Infrared Absolute Spectrometer (FIRAS) on-board the Cosmic 
Background Explorer (COBE). Shorter wavelengths have been investigated with the help 
of the COBE Diffuse Infrared Background Experiment (DIRBE) from 240 /im to 1.25 /xm. 
Hauser et al. (1998a) measured the CIB at 240 and 140 fim and set upper limit for other 
wavelengths. Lagache et al. (2000) proposed an estimation of the CIB at 100 /im. 

The CIB is the signal that remains after subtracting from the total celestial brightness 
the emission from the interstellar medium, stars and interplanetary dust (scattering and 
emission). However, estimating the contamination due to these three foreground components 
is challenging. The interstellar medium is the main contaminant for long wavelengths whereas 
starlight and zodiacal light dominate at wavelengths studied in the present paper (J — 
1.25 /xm and K = 2.2 //m). 

Stellar population models were used to estimate the brightness of starlight in the large 
beam DIRBE data (Arendt et al. 1998) but this technique led only to upper limits. De- 
tections of the CIB at 3.5 fim (L) and 2.2 fim (K) were reported by Gorjian et al. (2000) 
using near-infrared observations of a 2° x 2° dark spot to measure the brightness of stars in 
a few DIRBE beams. Wright (2001) confirmed the result at 2.2 //m and proposed a weak 
hmit for 1.25 /xm (J) using data from the 2 Micron All-Sky survey (2MASS) for 4 dark 
spots. Matsumoto et al. (2000) estimated the CIB from 1.4 /xm to 4 //m with the Near 
Infrared Spectrometer (NIRS) on-board the InfraRed Telescope in Space (IRTS). Although 
NIRS provided a point source catalog limited to stars brighter than 7.5 mag at 2.24 ^m, 
significant stellar contamination remains and the CIB estimate is strongly dependent on the 
stellar population model used to remove contribution from fainter stars. 

This work presents a large scale study of the CIB using 2MASS data for 1400 deg^ 

of the sky in order to estimate accurately the stellar contribution to the surface brightness 
observed by DIRBE. Section 2 describes the 2MASS and DIRBE data. The method used 
to compare these two data sets is presented is Section 3, and the following section presents 
the determination of the CIB for J and K with the associated uncertainties. Section 5 is 
dedicated to the comparison with previous results and with galaxy counts. 
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2. 



Data 



2.1. 



2MASS 



The 2MASS second Incremental Data Release from the 2MASS survey covers 48% of the 
sky and contains 1.6 x 10^ stars (Cutri et al. 2000). The Point Source Catalog completeness 
limits are 15.8 and 14.3 mag with a signal-to-noise ratio greater than 10:1 for J and Kg, 
respectively. No photometry is provided for stars brighter than 4-5 mag (see Sect. 3.2). 

Since the catalog of the whole release is ~ 7 GB in binary format, we choose to work 
on integrated maps in which each pixel corresponds to the integrated flux in a 5' x 5' box. 
These maps are constructed from the point source catalog clipped to remove sources fainter 
than the completeness limit. Magnitudes are converted to flux density using the flux for zero 
magnitude Fq{J) = 1603 Jy defined by Campins et al. (1985) and Fq{Ks) = 698 Jy which 
is an extrapolation of the J value to the Kg filter assuming a black body at 9700 K (Vega 
photometric system). Uncertainties on these photometric zero points are about 2% and do 
not affect our results. Each 5' x 5' box contains the sum of the point source fluxes in the 
area. 

The resulting image in an Aitoff (equal area) projection consists of 4321 x 2161 pixels. 
This size is reasonable and the image can be handled globally. Errors resulting from this 
pixelization are negligible since 5' is small compared to the DIRBE beam of 1° (see below). 

We should point out that the 2MASS Ks filter (2.17 fim) and the DIRBE K filter 
(2.2 yum) are slightly different. Persson et al. (1998) discuss the K to Kg transformation 
and present a list of photometric standards in both colors. Moreover, as Pahre & Mould 
(1994) point out, the Kg bandpass is intermediate between K and K' and rough corrections 
of Kg to K can be obtained by averaging between K' and K. This transformation is well 
defined by Wainscoat & Cowie (1992): K' - K ^ (0.22 ± 0.03) (i/ - i^T) for < i/ - i^T < 0.4 
(K'—K ~ 0.07 otherwise). Hence, for a K5 spectral type star we would obtain Kg — K ~ 0.01 
for example. This correction is negligible compared to uncertainties described further in the 
paper and differences between 2MASS and DIRBE filters will be ignored in the following 
analysis. 



Cryogenic DIRBE operation took place from 1989 December to 1990 September. During 
these 10 months, the sky was observed at the rate of half of the sky per week in 10 bands from 
1.25 to 240 fj,m (Hauser et al. 1998a). The DIRBE instrument was designed to make ac- 



2.2. 



DIRBE 
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curate absolute sky brightness measurements, with a strayhght rejection of < 1 nW m~^ sr~^ 
and an absolute brightness cahbration uncertainty of 0.05 and 0.03 nW m~^ sr~^, at 1.25 
and 2.2 fxm, respectively. 

The zodiacal light is the first component along the hght of sight and must be removed 
using a dust interplanetary model (e.g. Kelsall et al. 1998). However these models are not 
unique (see Wright 1998; Gorjian et al. 2000 for another model). Artifacts still remain. In 
the following analysis, we use the Zodi-Subtractcd Mission Average (ZSMA) maps produced 
by the DIRBE team (Hauser ct al. 1998b). The zodiacal light intensities were subtracted 
week by week and the residual intensity values were averaged to create the ZSMA maps. 

To compare DIRBE and 2MASS, we need a precise knowledge of the DIRBE beam in 
both J and K bands. The 42' beam size commonly mentioned in the literature is only valid 
for daily maps, not the annual maps used here. The beam for the J band is presented in the 
Figure 1 . This is an effective beam profile (provided with DIRBE maps) , which measures the 
relative response of DIRBE to a point source, including the effects of sky-scanning and data- 
sampling rates. The beam for the annual maps corresponds to an average of the effective 
beam for all orientations (Fig. 1). Similar results are obtained for J and K and the full 
width half maximum of both beams is ~ 1° (sohd angle ~ 0.78 deg^). 

3. Method 

3.1. 2MASS integration on DIRBE pixels 

First, we exclude from this study regions of the sky for which we know that the GIB can- 
not be straightforwardly extracted, for example due to high stellar density fields, interstellar 
cirrus (which is responsible for emission and scattering), and regions of residual zodiacal 
light structures in the DIRBE maps. 

We keep only regions that satisfy the following criteria: 

• high galactic latitude: |6| > 40° 

• low DIRBE 240 //m flux in order to eliminate diffuse emission or scattered fight from 
cirrus clouds: /240 < 3 MJy sr~^ 

• fiigfi ecliptic latitude: \(3\ > 30° 

• exclusion of the Magellanic clouds 
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The most straightforward way to remove the stellar component (2MASS) from the 
DIRBE flux is to work directly in the COBE coordinate system (projection on a cube) and 
to integrate the 2MASS maps on each DIRBE pixel. However, since the 2MASS data still 
contain coverage holes, we kept only those that are completely covered by 2MASS. The 
effective useful region covers finally about 1400 deg^ (Fig- 2). 



3.2. Bright stars 

In the remaining 1400 dcg^ piece of sky, ~ 25% of pixels are contaminated by bright 
stars. Stars bright enough to be detected as point sources in DIRBE are saturated in 
2MASS. Consequently their photometry cannot be derived from 2MASS to be subtracted 
from DIRBE. 

Moreover, the DIRBE beam has a size of 1°, but the sampling of the data is 22'. Several 
pixels are therefore contaminated by a bright star, and this contamination depends on the 
star position in the beam. To identify pixels affected by these stars we use a median filter 
which deals with the 2MASS holes. A pixel is considered as contaminated if it is brighter 
than the threshold, T, defined by the median plus 2(T, M + 2(t, of the surrounding pixels in 
a 3° radius circle of the DIRBE-2MASS map: 

T'^M'i DIRBE - 2MASS ) + SaM DIRBE - 2MASS ] (1) 

\«,r=3° / \i,r=3° / 

Pixels fainter than M — 2a are also removed to avoid bias. In order to validate this operation, 
we use the 2MASS Bright Infrared Star Compilation (Tam et al. 1999^) which contains 
photometry for stars as bright as —4 mag in JHKg. This compilation consists of data taken 
from the literature and photometry extrapolations from IRAS, MSX or 2MASS. Among the 
343 stars brighter than J — 5 mag in our regions of interest, 11 are not identified with 
the median filter (Eq.l). However, these remaining contaminated pixels are rejected by the 
robust hnear fitting of DIRBE/2MASS brightness (see Fig. 4). 

After removing bright stars, 1040 deg^ of the sky remain. 



^http://spider. ipac.caltech.edu/staff/raymond/2niass/birsc/birsc. html 
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Fig. 1.— Left: effective DIRBE beam for daily maps (FWHM= 42') at 1.25 yum. Right: 
averaged beam for annual map (FWHM= 1°). 




Fig. 2. — Aitoff projection in galactic coordinates of the integrated 2MASS map (J) on 
DIRBE pixels. The large hne represents the echptic plane. This study is hmited to the 
areas inside the dashed curves (i.e. high galactic and ecliptic latitudes) and non-zero pixels 
correspond to ~ 1400 deg^. 
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3.3. Faint stsirs 

The 2MASS integrated brightness maps were hmited to the formal completeness limits 
of 2MASS (j''™ = 15.8 mag, = 14.3 mag). A model is required to estimate the 

contribution of the fainter stellar population with respect to the position in the sky. We use 
a model in which the number and distribution of Milky Way stars as seen in the near to mid- 
infrared is adapted from the Bahcall & Soneira (1980) optical star count model. We have 
employed the discrete population formalism of Elias (1978), Jones et al. (1981), Garwood 
& Jones (1987) and Jarrett (1992). The model includes the class III (evolved giant), class 
IV (subdwarfs) and class V (main sequence) stellar populations. These are further divided 
into disk and spheroid spatial distributions (ala Bahcall & Soneira). The stars are discretely 
binned according to their spectral types, ranging from the hottest O stars to the coolest 
M dwarfs, giving a total of 22 separate spectral bins for the main sequence, and 12 for the 
luminous giants (ranging from G2 to M7 giant). The optical/infrared colors, fluxes and 
luminosity function per spectral type are based on empirical data (e.g., Koornneef 1983; 
2MASS, ISO and IRAS). Interstellar extinction is applied as a smooth exponential function 
of Galactic position, characterized by a scale height and disk length. The model parameters 
were tuned using deep optical and infrared star counts. The resultant model was validated 
using the 2MASS survey down to J = 15 and Kg — 14 mag (~ 2 mJy). 

The result of the model for a typical high galactic latitude field (/ = 210°, h — —50°) is 
presented in Figure 3. As expected, 2MASS counts are complete at least to Kg — 14.3 mag. 
According to the model, less than 3% (0.1%) of the total stellar energy is contained in stars 
fainter than Kg = 14.3 (19). The model is useful for magnitudes ranging from 14.3 to ~ 20 
in Kg. Assuming a conservative number density uncertainty of a factor of 2 for these faint 
magnitudes, the resulting error will not exceed 3% of the total stellar flux. 

4. Result 

4.1. DIRBE / 2MASS correlation 

We choose 4 areas with a radius of 5° in order to illustrate the correlation between 
DIRBE and 2MASS. Pixels contaminated by bright stars have been filtered. Figure 4 show 
the correlations for J and K bands. The main characteristics for each region are summarized 
in Table 1. We note that slopes are not exactly 1 (~ 1.07 for J, ~ 0.94 for K). This arises 
from the different calibration strategy used by DIRBE and 2MASS: DIRBE is calibrated with 
Sirius, 2MASS with a list of calibration stars of many different spectral types (Cutri et al. 
2000). These calibration differences affect only the slope of the DIRBE/2MASS relation, not 
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the value of the intercept which is the measure of the CIB. 

Correlation coefficients show that the correlation is good for both colors (see Table 1). 
The CIB is obtained by removing the faint star contribution to the intercept value. 

4.2. Zodiacal contamination 

Since the main error is suspected to come from the zodiacal subtraction in the DIRBE 
maps, a representation of the CIB versus the ecliptic latitude (3 is useful. For each pixel, we 
estimate the intercept of the DIRBE/2MASS linear correlation in a 5° radius area, and we 
remove the faint star model from the DIRBE data. After this operation we select regions 
with the best DIRBE/2MASS correlation as defined in Table 2. These values were obtained 
by examination of the histogram of each quantity, and correspond to the peak value divided 
by 4. This selection leads to ~ 550 deg^ of rehable regions in the sky. 

Figure 5 shows the resulting CIB versus ecliptic latitude for both colors with an ad- 
ditional plot of the DIRBE 25 yum Zodi-subtracted surface brightness. For this figure, we 
include both low and ecliptic latitude points. The 25 ^m band is the most sensitive to the 
zodiacal light and contains subtraction residuals from the interplanetary dust model. The 
comparison of the 25 ^m. plot with the CIB plots confirms that the scatter at low ecliptic 
latitude comes mainly from residual zodiacal effects. 

Systematic uncertainties associated with the zero point of the interplanetary dust models 
are difficult to estimate. Kclsall et al. (1998) estimate them by choosing the largest among 
three models at high ecliptic latitude. Accordingly, a conservative estimate of the systematic 
uncertainties in J and K is 6.25 kJy sr~^ and 4.4 kJy sr^^, respectively. If the three models 
used to derive these numbers were equally distributed in the space of possible models, the 
given uncertainties would correspond to 1.6(7. Unfortunately, the models are probably not 
equally probable and it is hard to interpret these numbers in terms of a confidence level. To 
be conservative, we assume they correspond to Icr. Theoretical uncertainties in the zodiacal 
fight model are further discussed by Dwek et al. (1998). 

Statistical uncertainties that comes from the DIRBE/2MASS correlation and from the 
zodiacal light can be estimated directly from Figure 5 by measuring the root mean square de- 
viation of the distribution. The resulting uncertainties are 2.7 kJy sr~^ for J and 2.1 kJy sr~^ 
for K. A summary of the different uncertainty contributions is presented in Table 3. The 
presence of fluctuations in the residual sky brightness, especially in J, suggests a small 
contribution to the sky brightness by starlight scattered off of interstellar dust grains. 
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5 10 15 20 

Ks 



Fig. 3. — Kg source counts from 2MASS (solid line), from the model used for faint stars 
(dashed line) and galaxy counts (diamond) from Gardner et al. (1997). The shaded area 
corresponds to the model uncertainties. The vertical line shows the cutoff used in the 2MASS 
catalog, the model provides us counts for fainter magnitudes. Note that at this location 
[l — 210°, h — —50°), galaxies dominate the counts for Kg > 16. 



Table 1. DIRBE/2MASS correlations for the four areas. 



band 


I 


h 


nb^ 


slope 


intercept^ 


corr'^ 


faint stars'* 


CIB^ 


J 


73.84 


70.02 


407 


1.00 ±0.03 


26.9 ± 1.6 


0.87 


1.44 


25.5 


K 






409 


0.89 ±0.02 


23.6 ± 1.1 


0.87 


1.84 


22.0 


J 


134.42 


50.30 


462 


1.08 ±0.03 


26.2 ± 1.7 


0.88 


1.68 


24.4 


K 






472 


0.93 ±0.02 


22.0 ± 1.0 


0.89 


2.34 


19.8 


J 


231.32 


-47.86 


386 


1.08 ±0.03 


21.2 ± 1.9 


0.88 


1.81 


19.2 


K 






389 


0.96 ±0.02 


20.3 ± 1.2 


0.91 


2.54 


17.9 


J 


319.54 


-84.42 


263 


0.98 ± 0.04 


26.8 ± 1.8 


0.82 


1.64 


25.2 


K 






260 


0.87 ±0.03 


23.9 ± 1.1 


0.86 


2.05 


22.1 



^number of pixels 

^expressed in kJy sr~^ 
■^correlation coefficient 

^faint stars (J > 15.8, Kg > 14.3) contribution from the Jarrett model in kJy sr 
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Fig. 4. — DIRBE versus 2MASS for 4 regions of 5° radius. Pixels contaminated by bright 
stars have been removed. Crosses correspond to pixels rejected by the Sex robust fitting. See 
Table 1 for main characteristics. 
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Fig. 5. — CIB versus ecliptic latitude for J (upper plot) and K (middle plot). Only pixels 
with galactic \h\ > 40°, good DIRBE/2MASS correlation and low DIRBE 240 ;um are plotted. 
The lower plot represents the total DIRBE flux at 25 /xm versus the ecliptic latitude for 
comparison. The same zodiacal light subtraction features appear in the 3 maps. The main 
uncertainty comes therefore from the zodiacal light and resulting CIB values are 22.9 ± 
7.0 kJy sr-i (54.0 ± 16.8 nW m'^ sr'^) and 20.4 ± 4.9 kJy sr^^ (27.8 ± 6.7 nW m^^ sr^i) 
for J and respectively. 
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4.3. 



Local interstellcir medium contamination 



Interstellar dust passing through the Solar System will scatter and reemit absorbed sun- 
light. The volume within which sunlight is significant (compared to the interstellar radiation 
field [ISRF]) is smaller than the size of the Solar System itself: at visible wavelengths, sunlight 
exceeds the ISRF out to the Oort cloud (~ 10'' AU), while at far-ultraviolet wavelengths, 
sunlight exceeds the ISRF only out to 800 AU (at 0.2 /xm) or less. The amount of interstellar 
dust close to the Sun was recently measured by the Ulysses and Gahleo spaceprobes (Griin 
et al. 1994); the volumetric cross-section is naCTd — 2.9 x 10^^'^ cm^ cm~^. It was found that 
the size distribution is deficient in small particles relative to the interstellar size distribution, 
such that the surface area is dominated by particles with size ~ 0.4 fj,m. 

To estimate the brightness of sunlight scattered by local interstellar dust, we assume 
the particles detected by Ulysses are spread uniformly throughout the Solar System. The 
emission from dust in the inner Solar System will have a detectable dependence on solar elon- 
gation angle, and would have been included in the DIRBE zodiacal emission model (Kelsall 
et al. 1998). Therefore, we integrate from 3 AU outward in calculating the brightness of sun- 
light scattered by interstellar dust and potentially contributing to an isotropic background. 
The brightness is 



where co is the albedo, Iq is the solar intensity, and $ is the scattering phase function. Using 
the properties of a mixture of astronomical silicates and graphite (Draine & Lee 1984), the 
albedo in the J and K bands is 0.42 and 0.21, respectively. For the phase function, we used 
a Henyey-Grccnstcin function with the asymmetry factor appropriate for the same mixture 
of astronomical silicates and graphite {gj = 0.15, gx = 0.02). The resulting brightness is 
-^lism — 0.8 and 0.3 kJy sr~^ in the J and K bands, respectively. 

Uncertainty in the estimate of the brightness of the local interstellar dust is due to three 
sources. First, our lack of accurate knowledge of the composition of the interstellar grains 
leads to a factor of approximately 2 uncertainty in their albedo; this estimate is based on the 
difference between the albedos of silicate and graphite grains (Draine & Lee 1984). Second, 
the lack of precise distinction between the emission that would be effectively isotropic and 
that which would be incorporated into the zodiacal light model leads to a ~ 20 % uncertainty; 
this is based on changing the minimum integration distance from 3 to 2.5 AU. Third, smaller 
particles that are not detected by Ulysses may be present at larger distances from the Sun; 
these particles are likely to be produce an insignificant scattering because of the weakness 
of sunlight at large distances from the Sun. 

To summarize, we estimate the GIB as follows: GIB( J) = 22.9 ± 7.0 kJy sr"^ (54.0 ± 
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16.8 nW m-2 sr-^) and CIB(K) = 20.4 ± 4.9 kJy sr^^ (27.8 ± 6.7 nW m-^ sr-^). 

5. Discussion 

5.1. Compctrison with previous studies 

Previous estimations of the CIB in K are 16.4 ±4.4 kJy sr"^ and 14.8 ±4.1 kJy sr~^ by 
Gorjian ct al. (2000) and Wright (2001), respectively. The difference between these values 
and the estimation given in this paper are within the measurement uncertainties (la) and 
thus are in good agreement. Most of the difference can be attributed to the zodiacal light 
subtraction because we used the Kelsall et al. (1998) model and they used the Wright (1998) 
model. 

The estimation of the J CIB can be compared to the weak limit of 12.0 ± 6.3 kJy sr~^ 
proposed by Wright (2001). The discrepancy with our value (22.9 ± 7.0 kJy sr~^) can be 
explained by the zodiacal light subtraction, since Wright would find 25.8 kJy sr~^ with the 
Kelsall et al. (1998) dust interplanetary model. 

Kashlinsky & Odenwald (2000) have investigated the CIB fluctuations in the DIRBE 
maps from 1.25 /xm to 4.9 /xm. Fluctuations on the DIRBE beam scale are expected to be 10- 
20% of the total CIB flux (Kashlinsky, private communication). Their values of and 
4.3127 kJy sr~^ (92% confidence level) for J and K respectively are therefore in agreement 
with our results in both bands. 



5.2. GalEtxy counts 

Galaxy counts and the CIB are closely related: at least part of the CIB is due to the 
integrated surface brightness of galaxies for all redshifts. In the past decade, many efforts 
have been made to obtain galaxy counts in a wide range of magnitudes. Jarrett & et al. 
(2001) performed large scale studies with 2MASS (2350 deg^) which provide galaxy counts 
for 8 < J < 15 and S < Kg < 14. Bershady et al. (1998) obtained ground based observations 
up to 24**^ magnitude in J and K with the Keck 10m telescope. Saracco et al. (1999) used 
ESQ NTT telescope and obtained galaxy counts for 17 < J < 24 and 16 < K < 22. Deeper 
observations are possible with the Hubble Deep Field (HDF) of the Hubble Space Telescope 
(HST) and Thompson et al. (1999) reached the 30*'^ magnitude in the HST filters F160W 
and FllOW which correspond roughly to / + J and H, respectively. 

Figure 6 shows the cumulative brightness of galaxies up to J = 30 and K — 30 mag. For 
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Table 2. Selection of the best DIRBE/2MASS correlation areas. 



Quantity 




J 




K 


number of pixels 


> 


80 


> 


80 


correlation coeff 


> 


0.72 


> 


0.82 


mean square dev. 


< 


200 


< 


80 


intercept a 


< 


2.8 


< 


1.6 


slope a 


< 


0.05 


< 


0.035 


slope 


G 


[0.97. 1.125] 


G 


[0.85.0.965] 



Table 3. Uncertainties (values for la). 



Quantity 



J 



K 



Bright Stars (< 5 mag) 
2MASS flux per DIRBE pixel 

(5 < J < 15.8, b<Ks< 14.3) 
Faint stars model 

(J > 15.8, Ks > 14.3) 
DIRBE flux 
DIRBE/2MASS slope 

CIB scatter 

(contains all stat. uncert.) 
Zodiacal light model 

(syst. uncert.) 
Local Interstellar Medium 
Total CIB uncertainty 



removed 

4.2 X 10"^ kJy sr" 

< 2 kJy sr"^ 
2.2 kJy sr-^ 



removed 

4.5 X 10"^ kJy sr" 

< 2 kJy sr"^ 
2.1 kJy sr"^ 



< 10%, no consequence < 10%, no consequence 
/ 2.5 X 10-^ kJy sr-i 



2.7 kJy sr-i 

6.25 kJy sr-i 
1.6 kJy sr~^ 
7.0 kJy sr-i 



2.1 kJy sr-i 

4.4 kJy sr~^ 
0.6 kJy sr~-^ 
4.9 kJy sr"^ 
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the faintest J and K region (magnitudes 24 to 30) we used the HST measurements assuming 
TTT^Fwow ~ = and niFiiow ~ '^k = 0. The resulting error can be neglected, since the 
contribution of these very faint galaxies to the total brightness is small (9% of the total 
energy in J, only 3% in K). But note that 90% of the energy is in the range 12 < J < 24 
and IA<K < 22. 
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Fig. 6. — Comparison of the CIB with the cumulative brightness derived from galaxy counts. 
Data are from 2MASS (Jarrett & et al. 2001), ESQ NTT (Saracco ct al. 1999), ESQ 2.2m 
(Saracco et al. 1997), Keck (Bershady et al. 1998) and HST (Thompson et al. 1999). Error 
bars correspond to a 95% {2a) confidence level. A significant part of the energy is still not 
detected in sources for both colors. 



Uncertainties on galaxy counts come from the Poisson statistics (cr = ^^N)^ and from 
completeness limitations and corrections. By comparing galaxy counts in the literature, we 
find that the number of galaxies by range of magnitude N is consistent to within a factor of 
1.75, which we adopt as a conservative value for the uncertainty. We find a J total brightness 
for galaxies of 3.6 ± 0.8 kJy sr^^, and a K total brightness of 5.3 ± 1.2 kJy sr~^, in good 
agreement with Pozzetti et al. (1998) who obtained 5.8 kJy sr~^ in the K band. 

The integrated galaxy counts are smaller than our DIRBE-2MASS CIB estimations, 
even compared with the 95% confidence level. We can address the possibility of large diffuse 
galaxy halos that would not be inchidcd in the galaxy luminosities (because of their faintness) 
and we find that to reduce the discrepancy between galaxy surface brightness and our CIB 
value to less than 1 a, an increase of each galaxy luminosity by a factor 3.6 and 2.4 is required 
in J and respectively. If we assume that galaxy photometries are correct, a significant 
part of the energy in the CIB is not detected in sources. The slope of the galaxy luminosity 
function is observed to flatten for K > 11 (Saracco et al. 1997) and the resulting slope 
for the cumulative brightness is almost for > 25 mag. For J, the slope change in the 
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luminosity function is not as pronounced asior K, but the slope of the cumulative brightness 
is also close to for J > 25. Consequently, the slope of galaxy luminosity function must 
increase for fainter magnitudes, or other contributions to the near-infrared background must 
be invoked. To summarize, the brightness of the cosmic background found in this paper 
(Sect. 4) is greater than the integrated brightness of galaxy counts by 19.3 ± 7.8 kJy sr~^ 
(46.3 ± 18.7 nW m'^ sr^^) at J and 15.1 ± 6.1 kJy sr^^ (20.6 ± 8.3 nW m'^ sr-^) at K. 

5.3. Other possible contributions to the CIB 

Although galaxies appear to be unable to explain the brightness of the cosmic back- 
ground, there are other possible sources of energy production that could lead to an isotropic 
cosmic background in the near-infrared. An early burst of star formation inside or outside 
of galaxies, (that is, either in primeval galaxies or in Population III stars), would result in 
a nearly uniform and isotropic background (Bond et al. 1986). The photospheric emission 
from these stars would be redshifted into the infrared. For ~ 10(100), the background 
due to high-mass stars will peak at wavelengths 1 (10) //m, and the background due to 
intermediate-mass stars will be in the range 0.5-2 (6) /xm. The gap between our background 
estimation and galaxy count contributions could be due to Population III stars. The present 
abundance of heavy elements limits the amount of early star formation such that Population 
III stars are not likely to be a very significant mass contribution to the universe as a whole 
(Bond et al. 1986), but the amount of star formation required to explain the near-infrared 
cosmic background does not violate metallicity constraints. The implications of such early 
star formation, for the small-scale density fluctuations of primordial matter and for the 
nature of the intergalactic medium, would be significant. 

Depending on the power spectrum of small-scale density fiuctuations, a population of 
black holes could form before galaxies. The accretion of surrounding material by these holes 
would generate energy that would make an isotropic background today. The temperature 
of the accretion disk would be ~ 10^ K, and the accretion luminosity would be redshifted 
into the near-infrared if the accretion is active at Zbh ~ 30-50 (Carr 1994). Our observed J 
band background brightness, after subtracting the integrated light from galaxies, could be 
explained by the accretion luminosity of a population of black holes. Such black holes would 
not comprise enough dark matter to close the universe, and their total mass would not violate 
constraints on the total baryon density inferred from models of primordial nucleosynthesis. 
But such a large population of objects (comparable to the total mass of galaxies) would 
contribute significantly to the total inventory of baryons. 

Massive particles that survive until redshift Zi, and subsequently decay into photons will 
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produce a background of decay photons that is subsequently redshifted (Bond et al. 1986). 

For example, part of the J and K band background could be produced by the decay of relic 
100 eV neutrinos at redshift Zi, ~ 300. The number of plausible models is very large, and we 
can only speculate whether the background is due to such decaying particles. 

6. Conclusion 

The integration of the 2MASS catalog over 550 deg^ (i.e. ~ 5500 DIRBE pixels) leads 
to an estimate of the CIB for J and K of 22.9 ± 7.0 kJy sr"^ (54.0 ± 16.8 nW m'^ sr"^) and 
20.4 ± 4.9 kJy sr"^ (27.8 ± 6.7 nW m'^ sr'^), respectively 

We have selected the most rehable areas of the sky by ehminating pixels contaminated 
by bright stars, low galactic latitude regions, low echptic latitude regions and possible cir- 
rus clouds contribution. A study of DIRBE— 2MASS versus ecliptic latitude clearly shows 
that the zodiacal light is responsible of most of the statistical uncertainty (~ 2 kJy sr^^). 
The interplanetary dust model contains a systematic error which is actually the dominant 
uncertainty in our results (~ 5 kJy sr~^). 

The integrated brightness of known galaxies (down to the SO*'^ magnitude) represents 
only a fraction of the CIB. Although the uncertainties due to residual contribution of the 
zodiacal light are significant, the discrepancy between DIRBE— 2MASS and galaxy counts 
would suggest that a significant part of the energy is still not observed in discrete sources. 
The slope of the galaxy luminosity functions may increase for galaxies fainter than the 30*^ 
magnitude, in which case the counts would agree with the CIB values. However, models 
predict other possible contributions to the background at these wavelengths (Bond et al. 
1986) such as a burst of star formation either in primeval galaxies or in Population III 
stars {z ~ 10), very massive black holes (accreting from an uniform pre-galactic medium at 
z 40), massive decaying big bang relic particles {z 300). We have put new constraints 
in the near-infrared that encourage revisiting the importance of those contributions to the 
CIB in cosmological models. 

The Space Infrared Telescope Facility (SIRTF) may extend this result for the L band 
at 3.5 jim. and might be able to constrain the contributions of the different components to 
the background. 

The final all-sky 2MASS catalogs will permit a more accurate study of the near-infrared 
excess: 100% of the sky will be still available after integration on the DIRBE beam, ~ 75% 
after bright star removal (only 18% for the 2MASS 2'^^'^ Release presently used) and about 
5000 deg^ will remain at high galactic and ecliptic latitudes. This coverage should allow us 
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to investigate structures in the zodiacal light and may help to improve interplanetary dust 
models. 
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